function y=loggamma(x,k,theta);

y=1/(theta^(k)*gamma(k))*exp(k*x-exp(x)/theta);
